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Abstract: This paper studies two stochastic BMAP arrival and bulk service C server queues (A) and (B) with k 
varying environments. The arrivals to the queue are governed by a batch Markovian arrival process of i version 
and the bulk service times are exponential with parameter U; in the environment i for 1 < i < k* respectively. 
When the environment changes from i to j, changes occur for arrival and service as follows: the arrival BMAP 
representation changes from the i version to the j version, the residual arrival time starts with the stationary 
probability vector of the j version BMAP, it becomes the initial j version upon arrival of customers and the 
exponential service time parameter changes from U; to Uj for 1 < i, j < k*. The system has infinite storing 
capacity and the service bulk sizes are finite valued random variables. Matrix partitioning method is used to 
study the models. In Model (A) the maximum of the arrival sizes M in all the environments is greater than the 
maximum of the service sizes N in all the environments, (M > N), and the infinitesimal generator is partitioned 
as blocks of the sum of the number of BMAP phases of all environments times the maximum of the arrival 
sizes for analysis. In Model (B) the maximum of the arrival sizes M in all the environments is less than the 
maximum of the service sizes N in all the environments, (M < N), where the infinitesimal generator is 
partitioned using blocks of the sum of the number of BMAP phases of all environment times the maximum of 
the service sizes for analysis. Five different cases associated with C, M and N due to partitions are treated. They 
are namely, (Al) M >N > C, (A2) M > C >N (A3) C >M >N, which come up in Model (A); (Bl) N > C and 
(B2) C >N, which come up in Model (B) respectively. For the cases when C < M or N Matrix Geometric results 
are obtained and for the cases when C > both M and N Modified Matrix Geometric results are presented. The 
basic system generator is seen as a block circulant matrix in all the cases. The stationary queue length 
probabilities, its expected values, its variances and probabilities of empty queue levels are derived for the 
models using Matrix Methods. Numerical examples are presented for illustration 

Keywords: Block Circulant, BMAP Arrival, Bulk Service, C servers, Infinitesimal Generator, Matrix methods. 

I. INTRODUCTION 

In this paper two multi server queues with batch Markovian arrival process (BMAP) and bulk service have been 
studied with random environments using matrix geometric methods. For M/M/l bulk queues with random 
environment models one may refer to Rama Ganesan, Ramshankar and Ramanarayanan [1] and M/M/C bulk 
queues with random environment models are of interest in Sandhya, Sundar, Rama, Ramshankar and 
Ramanarayanan [2]. PH/PH/1 bulk queues without variation of environments have been treated by Ramshankar, 
Rama Ganesan and Ramanarayanan [3] and the same type of queues with random variation of environments are 
studied by Ramshankar, Rama, Sandhya, Sundar and Ramanarayanan [4]. Bini, Latouche and Meini [5] have 
studied numerical methods for Markov chains. Chakravarthy and Neuts [6] have discussed in depth a multi- 
server queue model. Gaver, Jacobs and Latouche [7] have treated birth and death models with random 
environment. Latouche and Ramaswami [8] have analyzed Analytic methods. For matrix geometric methods 
and models one may refer Neuts [9]. Batch Markovian arrival processes are presented by Lucantony [10] and 
are analyzed also by Cordeiro and Khroufch [11]. The models considered in this paper are general compared to 
existing queue models in literature since a BMAP arrival and multi server bulk service queue with random 
environment has not been studied at any depth so far. The number of servers increases with the arrival of 
number of customers till it becomes C. Usually the partitions of the bulk arrival models have M/G/l upper - 
Heisenberg block matrix structure with zeros below the first sub diagonal. The decomposition of a Toeplitz sub 
matrix of the infinitesimal generator is required to find the stationary probability vector. In this paper the 
partitioning of the matrix is carried out in a way that the stationary probability vectors have a Matrix Geometric 
solution or a Modified Matrix Geometric solution for infinite capacity C server bulk arrival and bulk service 
queues with randomly varying environments. 
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Two models (A) and (B) on BMAP/M/C bulk queue systems under k* varying environments with infinite 
storage space for customers are studied here using the block partitioning method. The M/PH/1 and PH/M/C 
queues with random environments have been studied by Usha in [12] and [13] without bulk arrivals and bulk 
services. It has been noticed by Usha in [12, 13] that when the environment changes the remaining arrival and 
service times are to be completed in the new environment. The residual arrival time and the residual service time 
distributions in the new environment are to be considered at an arbitrary epoch since the spent arrival time and 
the spent service time have been in the previous environment with distinct sizes of PH phase. Further new 
arrival time and new service time from the start using initial PH distributions of the new environment cannot be 
considered since the arrival and the service have been partly completed in the previous environment indicating 
the stationary versions of the arrival and service distributions in the new environments are to be used for the 
completions of the residual arrival and service times in the new environment and on completion of the same the 
next arrival and service onwards they have initial versions of the PH distributions of the new environment. The 
stationary version of the distribution for residual time has been well explained in Qi-Ming He [14] where it is 
named as equilibrium PH distribution. Ramshankar, Rama, Sandhya, Sundar, Ramanarayanan in [4] have 
studied PH/PH/1 queue models with bulk arrival, bulk service with random environment introducing the 
stationary version for the residual times. In this paper the stationary probability starting vector of the new 
version is used when the environment changes for the residual arrival time and it becomes the initial new 
version of BMAP distribution after the arrival. Model (A) presents the case when M, the maximum of all the 
maximum arrival sizes in the environments is bigger than N, the maximum of all the maximum service sizes in 
all the environments. In Model (B), its dual, N is bigger than M, is treated. In general in Queue models, the state 
space of the system has the first co-ordinate indicating the number of customers in the system but here the 
customers in the system are grouped and considered as members of M sized blocks of customers (when M >N) 
or N sized blocks of customers (when N > M) for finding the rate matrix. For the C server system under 
consideration, Model (A) gives three cases namely (Al) M > N > C, (A2) M > C > N and (A3) C > M > N and 
Model (B) gives two cases namely (Bl) N > C, and (B2) C > N. The case M=N with various C values can be 
treated using Model (A) or Model (B). The matrices appearing as the basic system generators in these models 
due to block partitions are seen as block circulant matrices. The stationary probability of the number of 
customers waiting for service, the expected queue length, the variance and the probability of empty queue are 
derived for these models. Numerical cases are presented to illustrate their applications. The paper is organized in 
the following manner. In section II and section III the BMAP/M/C bulk service queues with randomly varying 
environment in which maximum arrival size M is greater than maximum service size N and the maximum 
arrival size M less than the maximum service size N are studied respectively with their sub cases. In section IV 
numerical cases are presented. 

II. MODEL (A). MAXIMUM ARRIVAL SIZE M GREATER 

THAN MAXIMUM SERVICE SIZE N 

2.1 Assumptions for M > N. 

(i) There are k* environments. The environment changes as per changes in a continuous time Markov chain with 
infinitesimal generator Q 1 of order k* with stationary probability vector §. 

(ii) In the environment i for 1 < i < k*, the batch arrivals occur in accordance with Batch Markovian Arrival 
Process with matrix representation for the rates of batches of size m < M ; given by the finite sequence {D m , 0 < 
m < } with phase order k t where D l 0 has negative diagonal elements and its other elements are non-negative; 
D l m has non-negative elements for 1 < m < M t . Let D l =Hm=o D m m d (p t be the stationary probability vector of 
the generator matrix D l with <pj) 1 = 0 and (p^ = 1 for 1 < i < k*. 

(iii) When the environment changes from i to j for 1 < i, j < k*, the arrival process BMAP of the j version starts 
as per stationary (equilibrium) probability vector of the j version of the arrival process for the completion of the 
residual arrival time there after the arrivals occur as per BMAP of the j version, namely, { D ] m 0 < m < M j }. 

(iv) Customers are served in batches of different bulk sizes. There are s servers to serve when s customers are 
present in the system for 1< s < C. When C or more than C customers are present in the system the number of 
servers to serve customers is C. In the environment i fori < i < k*, the time between consecutive bulk services 
has exponential distribution with parameter s/^ when s customers (s servers ) are in the system for 1< s < C and 
with parameter Cfi t when C or more than C customers (C servers )are present where [i t is the parameter of 
single server exponential service time disttibution. At each service epoch in the environment i, ip t customers are 
served with probabilities given by P (xpi = j) = q] for 1< j < JV; when more than N t customers are waiting for 
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service where Y lj . l =1 q l j =1. When n customers n < iV; are in the system, then j customers are served with 
probability, for 1< j < n-1 and n customers are served with probability^ l n q) for 1 < i < k*. 

(v) When the environment changes from i to j, the exponential service time parameter of single server changes 
from fa to fij , the bulk service size ip t changes to xpj and the maximum service size N t changes to Nj . 

(vi) The maximum batch arrival size of all BMAPs', M= max^^M; is greater than the maximum service size 
N= nra^^.JVj 

2.2.Analysis 

There are three sub cases for this model namely (Al) M > N > C, (A2) M > C >N and (A3) C > M >N. Sub 
Cases (Al) and (A2) admit Matrix Geometric solutions and they are treated in sub section (2.2.1). Modified 
Matrix Geometric solution is presented for Sub Case (A3) which is studied in sub section (2.2.2). The state of 
the system of the continuous time Markov chain X (t) under consideration is presented as follows. 
X (t) = {(n, m, i, j): for 0 < m < M-l; 1 < i < k*, 1 < j < k l and n > 0}(1) 
The chain is in the state (n, m, i, j) when the number of customers in the system is n M + m, for 0 < m < M-l, 0 
< n < co, the environment is i for 1 < i < k* and the arrival phase is j for 1 < j < k t . When the number of 
customers in the system is r, then r is identified with (n, m) where r on division by M gives n as the quotient and 
m as the remainder. . Let the survivor probabilities of services xp t be respectively for the environment state i for 
1 < i < k*. P{rp i >m)= Q l m =\-^n=i<ln, for 1 < m < iV ; -1 (2) 
Q l m =0 for m > JV ; and Q l 0 = 1. (3) 

2.2.1 Sub Cases: (Al) M > N > C and (A2) M > C > N 

When M>N>CorM>C>N, the BMAP/M/C bulk queue admits matrix geometric solution as follows. The 
chain X (t) describing them, has the infinitesimal generator Q A21 of infinite order which can be presented in 
block partitioned form given below. 

" A 0 0 0 . . . 

A, An 0 . . 



0 0 



0 



In (4) the states of the matrices are listed lexicographically as 0, 1, 2, 3, .... n, .... Here the vector n is of type 1 x 
ME^iti and n = ( (n, 0, 1, l),(n, 0, 1, 2)...(n,0,l, k t \ (n, 0, 2, l),(n, 0, 2, 2)...(n,0,2, k 2 ),. . .,(n, 0, k*, 
l),(n, 0, k*, 2)...(n, 0, k*, k kt ), (n, 1, 1, l),(n, 1, 1, 2)...(n, 1, 1, k t ), (n, 1, 2, l),(n, 1, 2, 2)...(n, 1, 2, 
fe 2 ),...,(n, 1, k*, l),(n, 1, k*, 2)...(n, 1, k*, k kt ),..., (n, M-l, 1, l),(n, M-l, 1, 2)...(n, M-l, 1, /q), (n, M-l, 2, 
l),(n, M-l, 2, 2)...(n, M-l, 2, fe 2 ),...,(n, M-l, k*, l),(n, M-l, k*, 2)...(n, M-l, k*, k kt ) ) for n > 0. The 
matrices B 1 and A x have negative diagonal elements, they are of order M£f=! k t and their off diagonal elements 
are non- negative. The matrices A 0 ,andA 2 have nonnegative elements and are of order M£f=! k t and they are 
given below. 

Let Q[= D l 0 + {-Ctii for 1 < i < k* (5) 

where I indicates the identity matrix of order given in the suffix, Q l is of order k t . Considering the change of 
environment switches on stationary version of BMAP arrival in the new environment, the following matrix Q of 
order £f=i k t is defined which is concerned with change of environment during arrival time and service time. 

r q'i 



•02,1 Q'2 



n 3 . 



(6) 



ln ktil n kti2 n kti3 - Q' k J 

where H^j is a rectangular matrix of type k t x kj whose all rows are equal to (Qi)jj <Pj for i ^ j , 1 < i, j < k*. In 
the environment i, for 1 < i < k*, the matrix of arrival rates of n customers corresponding to the arrival in 
BMAP is D l n which is a matrix with non-negative elements for 1 < n < M t and D l n = 0 matrix for n > M t (7) 
and the rate with which n customers are served by a single server for 1< n < iV; is given by 
Sln^iln and s 'i,n= 0 if n > iV;. (8) 
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Let A n 



for ! 



: M 



(9) 



Lo 0 0 ••• D**J 

In (9) A n is a square matrix of order YJt=t k t ; D^is a square matrix of order kj for 1 < j < k* and 0 appearing as 

(i,j) component of (9) is a block zero rectangular matrix of type fcjXfc,. 

S' ln I kl 0 0 ••• 0 



Let U„ ■■ 



S 2 ,nh 2 
0 



0 



•3,n'k 3 



for 1 <n<N 



(10) 



» o 0 s ktin i kk 

In (10) U n is a square matrix of order£f=ifc;; Sj n I k . is a square matrix of order fc ; for 1 < j < k* and 0 
appearing as (i, j) component of (10) is a block zero rectangular matrix of type k t x kj . The matrix A t for i = 0, 
1, 2 are as follows. 



r a m 



o 



o 



0 



0 1 



0 CU N CU N _ 1 
0 0 CU N 



CU N 





r n 


Ai 


^2 




Am-N-2 


A M -n- 


1 ^M-W 


- A M _ 2 








n 






Am-N-2 


A M N- 


2 A M N i 


- A M _ 3 






cu 2 


CUi 


n 




Am-N-4 


A M N- 


3 ^M- W -2 


••• A M _ 4 


A M - 3 




CU N 




CJ/jV-2 




n 


Ai 




••• A M _ N _ 2 


A M -N-1 


A x = 


0 


CU N 






cu x 


n 




••• Am-n-3 


A M -N-2 




0 


0 






cu 2 


CDi 


n 


••• ^M-iV-4 


A M -N-3 




0 


0 


0 




CU N 




1 cf/ w _ 2 


n 


A t 




. 0 


0 


0 




0 








a 


For 


defining 


the matrices 


Bi 


the 


following 


component 


matrices 


Using (2) and (3) let 


7 'i,n = 




<n<N 


-1 which is a matrix of order k t for 1 ; 




V' ln 


0 


0 


0 












% = 


0 


v'2,,1 


0 ■•• 


0 


for 1 <n<N-l. 










0 


0 


0 - 


V'k*,n 













required 



,14, 

This is a matrix of order (£?=i fc;) and 0 appearing in the (i, j) component is a 0 matrix of type k t x fe ; for 1 < i, 
j<k*. 

n ... n 

(15) 



LetU = 



0 
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In (15), U is matrix of order (Hf=i fe;) and 0 appearing in the (i, j) component is a rectangular 0 matrix of type 
fcj x fey for 1 < i, j < k*. To write ^the block for 0 is to be considered which has queue length L— 0, 1, 2. . .M-l. 
When L = 0 there is only arrival process without service. The change in the environment from i to j switches on 
BMAP j version as per stationary (equilibrium) probability vector in the new environment j whenever it occurs 
for 1 < i, j, < k*. In the empty queue (L=0) when an arrival occurs in the environment i both the arrival time and 
the service time start. In block 0 when L =1, 2,..., M-l all the processes arrival, service and environment are 
active as in other blocks n for n > 0. Considering the change of environment switches on BMAP arrival process 
in the new environment through the stationary (equilibrium) probability vector when the queue is empty, the 
following matrix O' of order £?=i fe;is defined which is concerned with the change of environment during 
arrival time and is similar to O defined in (6). 

■ t\ n 12 n 13 ••■ n lkt 



n 3 , 



T' 3 



(16) 



ln ktil n k%2 n kti3 - r k J 

Here T\= D l 0 + diag(Q 1 ) i i and i2 i; - is a rectangular matrix of type fc ; x fe ; whose all rows are equal to (Qi)ij 
(Pj presenting the rates of changing to phases in the new environment for i ^ j and 1 < i, j < k*. 
The matrix B 1 for Sub Case (Al) where N > C and Sub Case (A2) where C > N are given below in (17) and 
(18) respectively. For the case when C=N, the matrix^ may be written by placing C in place of N in the N-th 
block row in (18) and there after the multiplier of Uj is C. Let Q x j = D.' - jU for 0 < j < C and Q 1C =Q 



CV c .j C£/ c _, CV c - 2 CU c - 3 

CVr CUr Ci/V-! CI/;-; 



;:• _ : :v _- 



CU:--:-- CV.-- : CUy-:-: 



(■V + l)£/„-, GV+1)J/ K - 
GV + 2)£/ K {X+1)Uk- 



For the case when M = C, the multiplier C does not appeal" as a multiplier for the Uj matrices in the matrix B- in 
(18) in the 0 block of (4) and C appears as a multiplier for all Uj matrices in the matrices of A- and A 2 from row 
block 1 onwards. The basic generator of the bulk queue which is concerned with only the arrival and service is 
a matrix of order [M£f=ife; ] given below in (21) where Q A =A 0 + A- + A 2 (19) 
Its probability vector w' gives, w'Q A =0 and w'. e = 1 (20) 
It is well known that a square matrix in which each row (after the first) has the elements of the previous row 
shifted cyclically one place right, is called a circulant matrix. It is very interesting to note that the matrix Q A is 
a block circulant matrix where each block matrix is rotated one block to the right relative to the preceding block 
partition. 
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In (21), the first block-row of type [£?=ifc t ] x[Mg=iM is, W = (£1 + A M ,A X , A 2 , 
■ ■■,A M _ N _ 2 , A M _ N _ 1 , A M _ N + CU N , A M _ 2 + CU 2 , A M _ X + CUi) which gives as the sum of the blocks 
(n+A M ) + A 1 +A 2 +...+A M _ N _ 2 + A M _ N _ X +A M _ N + CU N +...+A M _ 2 + CU 2 + A M _ X + CU t = O" which is 
the matrix given by 



Q"i 
n 2 ,i Q"; 



(22! 



where using (5) and (6), Q'\ = D l + diag (Q 1 ) £ £ for 1 < i < k*. The stationary probability vector of the basic 
generator given in (21) is required to get the stability condition. Consider the vector 
w = (4>i<Pi$2 ( P2,--,<bk*<Pk*) (23) 
where (|> = (4> 1( 4> 2 , ■■■ , 4>fc*) is the stationary probability vector of the environment, (p t = ( <p i; ) is the stationary 
probability vector of the arrival BMAP D l for 1 < i < k*. It may be noted (p^fi 1 ^. This gives <p i (p i Q" i = 
$;C2i);,; <Pi for 1 < i < k*. Now the first column of the matrix multiplication of wQ" is $i (<2i)i,i<Pi,i + 

< Pii.Q\)2,\ < P\\\. < P2 e '\ + + < Pk*(Qi)k*,i ( Pn l<Pk* e ] = 0 since (<Pj)e = 1 and ^>Q X =0. In a similar manner it can be 

seen that the first column block of the matrix multiplication of wQ" is 4*i(Qi)i,i^i + 

4>2 (Qi)2,i ( Pi [( < P2) e l + + 4>/c* (Qi)k*,i ( Pi[( ( Pk*) e ] = 0 and i-th column block is 

<\>i{Qi)u ( Pi[{<Pi)e}+<\> 2 {Q 1 ) u <p l [{<p 2 )e}+ +<b i (.Q 1 ) i , t <p i +...+<b k ,(Q 1 ) k *, i <p i K<p k M= 0. This shows that 

w(f2 +A M ) + wA t + wA 2 +...+wA M _ N _ 2 + wA M _ N _ 1 + wA M _ N + wCU N +...+wA M _ 2 + wCU 2 + wA M _ 1 + 
wCt/ 1 = w Q"=0. So (w, w,...,w) .W= 0 = (w, w, ....w) W where W is the transpose W. This shows 
(w, w...w) is the left eigen vector of Q' A and the corresponding probability vector is 

= h^L^L HL\ (24) 

\M MM MJ v ' 

where w is given by (23). Neuts [5], gives the stability condition as, w' A 0 e < w' A 2 e where w' is given by 
(24). Taking the sum cross diagonally in the A 0 and A 2 matrices, it can be seen using (9), (10), (11) and (12) 
that 

w' A 0 e=± WQXU nA n )e=± (£^ =1 Zf=i n * ( (p t ^)e ) = 0, Sf =1 n (^)e 

^(ZU^tCPi&^nDtie <w'A 2 e=± w(^ =1 nCUje=^ (l^I^n^/i^e) ) = 
jj (Sf=i (pilJ-i Sn=i nq l n )) =^(2*=i <PiHiE(ipi) . This gives the stability condition as 

I'[*i<l>i<Pi(Z%=inD l n )e<Cp2 1 (25) 
This is the stability condition for the BMAP/M/C bulk service queue with random environment for 
Sub Case (Al) M > N > C and Sub Case (A2) M > C > N. When (25) is satisfied, the stationary distribution 
exists as proved in Neuts [9]. Let ji (n, m, i, j), for 0 < m < M-l, 1 < i < k*, 1 < j < k t and 0 < n < oo be the 
stationary probability of the states in (1) and 7T n be the vector of type lxM£f=! k t with n n = (?i(n, 0, 1, 1), 7t(n, 
0, 1, 2) ... 7i(n, 0, 1, k x ), 7i(n, 0, 2, 1), n(n, 0, 2, 2),...7t(n, 0,2, k 2 )... 7t(n, 0, k*, 1), 7i(n, 0, k*, 2),...n(n, 0,k*, 

k kt ) 7t(n, M-l, 1, 1), 7i(n, M-l, 1, 2) ... 7t(n, M-l, 1, /q), 7i(n, M-l, 2, 1), 7t(n, M-l, 2, 2),...7i(n, M-l, 2, 

k 2 )... 7i(n, M-l, k*, 1), 3t(n, M-l, k*, 2),...jt(n, M-l,k*, k kt )for n > 0. The stationary probability vector n = 

(77:0,77^,713, ) satisfies nQ A 21 =0 and 7re=l. (26) 

From (26), it can be seen n 0 B 1 + 71^=0. (27) 
Ti n _ 1 A 0 +n n A 1 +n n+1 A 2 = 0, for n > 1. (28) 
Introducing the rate matrix R as the minimal non-negative solution of the non-linear matrix equation 
Aq+RA^R 2 A 2 =0, (29) 
it can be proved (Neuts [9]) that n n satisfies 7T n = 7T 0 R n for n > 1. (30) 

Using (27) and (30), tt 0 satisfies tt 0 [B t + RA 2 ] =0 (31) 

The vector 7T 0 can be calculated up to multiplicative constant by (31). From (26) and (30) 
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7r 0 (/-/?)- 1 e=l. (32) 
Replacing the first column of the matrix multiplier of n 0 in equation (31) by the column vector multiplier of n 0 
in (32), a matrix which is invertible may be obtained. The first row of the inverse of that same matrix is n 0 and 
this gives along with (30) all the stationary probabilities. The matrix R given in (29) is computed using 
recurrence relation R(0) = 0; R(n + 1) = -A^A^ 1 -R 2 (n)A 2 Aj 1 , n > 0. (33) 
The iteration may be terminated to get a solution of R at an approximate level where ||i?(n + 1) — /?(n)|| < s. 
2.2.2 Sub Case: (A3) C > M > N 

When C > M > N, the BMAP/M/C bulk queue admits a modified matrix geometric solution as follows. The 
chain X (t) describing this Sub Case (A3), can be defined as in (1) presented for Sub Cases (Al) and (A2). It has 
the infinitesimal generator Q A2 2 °f infinite order which can be presented in block partitioned form given below. 
When C > M, let C = m* M + n* where m* is positive integer and n* is nonnegative integer with 0 < n* < M-l. 
[B\ A 0 0 0 0 ••• 0 0 0 0 
An A 0 0 0 ••• 0 0 0 0 




In (34) the states of the matrices are listed lexicographically as 0, 1, 2, 3, .... n, .... Here the vector n is of type 
lxMEj^fcjand n = ( (n, 0, 1, l),(n, 0, 1, 2)...(n, 0, 1, /q), (n, 0, 2, l),(n, 0, 2, 2)...(n, 0, 2, fc 2 ),...,(n, 0, k*, 
l),(n, 0, k*, 2)...(n, 0, k*, k kt ), (n, 1, 1, l),(n, 1, 1, 2)...(n, 1, 1, /q), (n, 1, 2, l),(n, 1, 2, 2)...(n, 1, 2, 
fe 2 ),...,(n, 1, k*, l),(n, 1, k*, 2)...(n, 1, k*, k kt ),..., (n, M-l, 1, l),(n, M-l, 1, 2)...(n, M-l, 1, /q), (n, M-l, 2, 
l),(n, M-l, 2, 2)...(n, M-l, 2, fc 2 ),...,(n, M-l, k*, l),(n, M-l, k*, 2)...(n, M-l, k*, k kt ) ) for n > O.The matrices 
B\,A X j fori < j < m * and A x have negative diagonal elements, they are of order M^fjq fc ; and their off 
diagonal elements are non- negative. The matrices A 0 andA 2 have nonnegative elements and are of order 
M£?=i fc; and the matrices A 0 ,A 1 and A 2 are same as defined earlier for Sub Cases (Al) and (A2) in equations 
(11), (12) and (13). Since C > M the number of servers in the system s equals the number of customers in the 
system L up to customer length becomes C. Once number of customers becomes L > C, the number of servers in 
the system remains C. When the number of customers becomes less than C, the number of servers reduces and 
equals the number of customers. The matrix A 2 j for 1 < j < m*-l is given below. 

" ■•• 0 )MU N jMUun jMU 2 " jMU 1 

••• 0 0 (jM + l)U N ••• 0'M+l)i/ 3 (JM + 1)U 2 



0 



0 



0 



0 



0 



0 



(JM + N - 2)U N (JM + N - 2)U N n 
0 (JM + N - 1)U N 



0 



0 



(35s 



The matrix A 2 mt is as follows given in (36) when C = m*M + n* and n* is such that 0 < n* < N-l. Here the 
multiplier of Uj in the row block increases by one till the multiplier becomes C = m*M + n* and there after the 
multiplier is C for Uj for all blocks. When N < n* < M-l, A 2 mt is same as in (35) for j = m* 
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33 jr.?r?.-io-j2 rae 5:. 15:. C; •= r.'-;l":'c 

■ ft.* * 

(jM + 2)U 2 (jM + 2) ffi.;«-2 
QM + X)U S (/Jf+JV)P*-i (/Jf+JV)P«-i 



M-2)£/.v (Af-2)U. V _ : 

o (M-i)y. v 



Q'lM-2 *t 

• /. . i- v folio-" 



0M+.V+ l)Ui 



ei,«-.v- : 



0 0 0 (jM + M-2)U K (jM + M-Z)U K -i - ei.,M- M -: -1 

[ o o o—o (jm+m-i)u n • UM+ii-iWi e liH 

The matrix A lmt is as follows when C = m*M + n* and n* is such that 0 < n* < N-l. The multiplier of Uj 
increases by one till it becomes C = m*M + n* and thereafter in all the blocks the multiplier of Uj is C. 



(Mm « +l)U l 




A t 


•• A K-M-1 A M-M-7 






(Mm - +2)£/ 2 


(Jm - +2)^ 




" A M - H - 4 Am-s-1 
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(Mm - +n« -i)U n .-i 


(jm - +n « -l)y R „ 






















CU S 




cu K _ 2 


■■ Sic i 




A M .„. 


0 


cu s 




- q[ c 


A H-.<C-3 


A M-K- 


0 


c 


0 


■■ cu„ cu..., 


•■ «ic 


- 1 : 


0 


0 


0 


0 CU H 




Cm 



.55! 



When n* = N or n* > N then, in the matrix A lmt , there is slight change in the elements. When n* = N, in the 
N+l block row and thereafter C appears as multiplier of Uj , and when n* > N with n* = N + rfor 1 <r< M-N-l, 
in the n*+l block row U N appears in the r + 1 column block. C appears as multiplier for it and as the multiplier 
of Uj thereafter in all row blocks respectively. The basic system generator for this Sub Case is same as (21) with 
probability vector as given in (24). The stability condition is as presented in (25). Once the stability condition is 
satisfied the stationary probability vector exists by Neuts [9]. As in the previous Sub Cases, 
^Qa,2.2=° and ne=1 - (4°) 

The following may be noted. Ti n ^o +n n+i^i +n n+2^2 = 0, for n > m*, the rate matrix R is same as in previous 
Sub Cases with same iterative method for solving the same and n n satisfies n n = n mt R n ~ m * for n > m*. (41) 
The set of equations available from (40) are 7r 0 B' 1 +7r 1 v4 21 = 0, (42) 

n i AQ+n i+l A li+1 +n i+2 A 2i i + 2 0, for 0 < i < m*-2 (43) 

and n m*-l^ L 0 +ll m* A -l,m* +n m*+l A 2 = 0. (44) 

The equation 7re=l in (40) gives £™*o~ l7r ; e + n m*(J-R)~ le = 1 (45) 
Using n mt+1 =n mt R and equations (42), (43), (44) and (45) the following matrix equations can be seen where 
Q'a,2.2 is g iven b y (48). 

(TTq,^,^, TT mt )Q' A,2.2=° (46) 

(rc 0 .%.% n mt ) [(/ _ jR )-i e ]=l (47) 
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0 



0 



0 



0 



- 0 0 0 0 0 ••• A 2imt RA 2 +A lm 

Equations (46) and (47) may be used for finding (n 0 ,n 1 ,n 3 , Tt m *)- Replacing the first column of the first 

column- block in the matrix given by (48) by the column vector multiplier in (47) a matrix which is invertible 

can be obtained. The first row of the inverse matrix gives (n 0 ,n 1 ,n 3 , 7T mt ).This together with equation (41) 

give all the probability vectors for this Sub Case. 
2.3. Performance Measures 

(1) The probability P(L = r), of the queue length L = r, can be seen as follows. Let n > 0 and j for 0 < m < M-l 
be non-negative integers such that r = n M + m. Then it is noted that 
P (L=r) =£t=i Syii n (n, m, i,j ), where r = M n + m. 

(2) P (Queue length is 0) = P (L=0) = E?=i n (0, 0, i,j ). 

(3) The expected queue level E(L), can be calculated as follows. 
For Sub Cases (Al) and (A2) it may be seen as follows. Since n (n, m,i,])=P[L = Mn + m, and environment 
state = i, arrival BMAP phase=j], for n > 0, 0 < m < M-l, 1 < j < k t and 1 < i < k*, 
E(L) =£™ =0 2m=o H^itfUnin.m.iJ) (Mn + m) =Yn=o^n- (Mn... Mn, Mn+1... Mn+1, Mn+2. . .Mn+2. . . 
Mn+M-1... Mn+M-1) where in the multiplier vector Mn appears £?=i k t times; Mn+1 appears £?=i k t times; 
and so on and finally Mn+M- 1 appears £?=i fc; times. So E(L) =M£^ =0 nn n e +n 0 (I — i?) _1 f . Here § is a 
(M^^kj ) xl type column vector £,=(0, ...0,1, ... ,1,2, ... ,2, ... ,M - 1, ...,M - 1)' where 0,1, 2,... M-l appear 
Y, k i=i k i times in order.This gives E (L) = n 0 (I - Ry 1 ^ + Mn 0 (I - R)~ 2 Re . (49) 
For Sub Case (A3), E(L) =l£ =0 Y.m=o £?=i Z* l =1 n (n, m, i,j ) (Mn + m) = MZ™ =0 nn n e + Y.n=o \f = 
MSn=o nn n e+X™ o _1 71 A + 7r m*( i "R)~ 1 ^- Letting the generating function of probability vector O(s) = T.T=o n i s ^ 
it can be seen, ®(s) = Yn= *o _1 t^s' +n m * s m *(I-Rs) _1 and £n=o nn n e = <D'(l)e = Yn= *o _1 ™; e+7r mt m*(I-R) _1 e 
+ 7T mt (I-R)~ 2 Re. Using this, it is noted that 
E(L) = M £™*(f 1 in t e + 7r mt m*(I-R)- 1 e + 7r mt (I-R)- 2 Re] o" 1 n& + ^(FR)" 1 ^ (50) 

(4) Variance of queue level can be seen using Var (L) = E (L 2 ) - E(L) 2 . Let n be column vector 
iH0,..0,l 2 ,...l 2 2 2 ,..2 2 ,... (M-l) 2 ,... (M-l) 2 ]' of type (ME?=i k t ) xl where 0,1, 2,... M-l appear 
5^=1 k t times in order. Then it can be seen that the second moment, for Sub Cases (Al) and (A2) 
E (Z 2 ) = En=oEm=o £?=i£yi 1 Jr(n,77i,i,;) [Mn + m] 2 =M 2 [l£ =1 n(n - l)n n e + Y.n=o nn n e] + 

Y.n=Q^nV + 2MX™ =0 n7T n f. 

So, E(L 2 )=M 2 [tt 0 (/ - R)~ 3 2R 2 e + n 0 (I - Ry 2 Re]+n 0 (I - i?)" 1 ^ + 2M n 0 (I - Ry 2 R% (51) 
Using (49) and (51) the variance can be written for Sub Cases (Al) and (A2). 
For the Sub Case (A3) the second moment can be seen as follows. 

E (Z 2 ) = l£=o 1Z~ = 1 YjU n (n, m, i,j ) [Mn + m] 2 = M 2 E£ =1 n(n - l)n n e + £? =0 ™i n e] + 
I^ =0 7r„77 + 2MS^ =0 nir n f = M 2 [<D"(l)e +X™o~ 1 in t e+7r mt m*(I-R)- 1 e + 7r mt (I-R)- 2 Re] + S^f 1 ^ + 
^(FR)" 1 ^ + 2M [ Y.T= *o _1 i^i f+7r mt m*(I-R)-^ + 7r mt (I-R)- 2 R §]. This gives 

E (L 2 ) = M 2 [ S™*! -1 i(t - 1)71; e + m*(m*-l)7T mt (/ — i?) _1 e +2m*7T mt (I-R)" 2 Re +27T mt (I-R)- 3 R 2 e 
+I™o _1 i^i e+7r m ,m*(I-R)- 1 e + 7r mt (I-R)- 2 Re] + T,™^ 1 + ^(FR)" 1 ^ +2M [ I^f 1 itt; f+7r mt m*(I- 
R)" 1 ^ + 7r mt (I-R)- 2 R ^. (52) 
(52) Using (50) and (52) the variance can be written for Sub Case (A3). 



III. MODEL (B). MAXIMUM ARRIVAL SIZE M LESS THAN 

MAXIMUM SERVICE SIZE N 

In this Model (B) the dual case of Model (A), namely the case, M < N is treated. Here the partitioning 
matrices are of order N£f=i k t and the customers are considered as members of N blocks. M plays no role in the 
partition where as it played the major role in Model (A). Two Sub Cases namely (Bl) N > C and (B2) C > N 
come up in the Model (B). (When M =N and for various values of C greater than them, or less than them or 
equal to them, both Models (A) and (B) are applicable and one can use any one of them.) The assumption (vi) of 
Model (A) is modified without changing others. 
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3.1 Assumption. 

(vi) The maximum batch arrival size of all BMAPs', M= maar 1 ^< fcj ,M i is greater than the maximum service size 

N=max 1 < i ^ ct iV i . 

3.2.Analysis 

Since this model is dual, the analysis is similar to that of Model (A). The differences are noted below. The state 
space of the chain is as follows defined in a similar way presented for Model (A). 
X (t) = {(n, m, i, j): for 0 < m < N-l, for 1 < i < k* for 1 < j < fc f and 0 < n < oo}. (53) 
The chain is in the state (n, m, i, j) when the number of customers in the queue is, n N + m, the environment 
state is i and the BMAP arrival phase is j for 0 < m < N-l, for 1 < i < k*,for 1 < j < k t and 0 < n < oo. When the 
customers in the system is r then r is identified with (n, m) where r on division by N gives n as the quotient and 
m as the remainder. 
3.2.1 Sub Case: (Bl) N > C 

The infinitesimal generator Q B3 ^ of the Sub Case (Bl) of Model (B) has the same block partitioned structure 
given in (4) for the Sub Cases (Al) and (A2) of Model (A) but the inner matrices are of different orders and 
elements. 

0 
0 



A\ A\ 
A" 2 A'\ A'\ 



0 



0 



In (54) the states of the matrices are listed lexicographically as 0, 1, 2, 3, ... . n, .... Here the vector n is of type 
lxNl^ii andn = ((n, 0, 1, l),(n, 0, 1, 2)...(n, 0, 1, k x ), (n, 0, 2, l),(n, 0, 2, 2)...(n, 0, 2, k 2 ),. . .,(n, 0, k*, 
l),(n, 0, k*, 2)...(n, 0, k*, k kt ), (n, 1, 1, l),(n, 1, 1, 2)...(n, 1, 1, k t ), (n, 1, 2, l),(n, 1, 2, 2)...(n, 1, 2, 
fc 2 ),...,(n, 1, k*, l),(n, 1, k*, 2)...(n, 1, k*, k kt ),..., (n, N-l, 1, l),(n, N-l, 1, 2)...(n, N-l, 1, /q), (n, N-l, 2, 
l),(n,N-l,2, 2)...(n,N-l, 2, fe 2 ),...,(n, N-l,k*, l),(n, N-l, k* 2)...(n, N-l, k* k kt ) ) for n > 0. 
The matrices, B'\, A" 0 ,A'\ and A" 2 are all of order N 2f=i ^£ - The matrices B'\ andA" 1 have negative 
diagonal elements and their off diagonal elements are non- negative. The matrices A" 0 and A" 2 have 
nonnegative elements. They are all given below. Using the same matrices presented in model (A), for Q, 
, U, Q' and Q 1J in (6), (9), (10), (14) to (17) the partitioning matrices are defined below. 



0 



CU„ CU n _-l 



0 0 



cu 3 
cu^ 
cu 5 



cu 4 



cu N _ M _ l cu N _ M _ 2 

CU N _ M CUm.! 

cu N _ M+1 CU N _ M 



CU N _- 
CU N _ 



cu N _ 3 
cu N _ 2 



cu N _ M _ 2 

CUn-m-i 
CU N _, 

cu N _ 3 



cu N _ M _ : 



cu N _ M _ 3 



CU N _, 
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C" - --. 




n A t 








CV'v-3 


cu K _ 3 


CI' -. 


- cu N _„_ 2 


cu v _ u _ 3 cu K . u . z 




Q 


.1 


• CVs-2 


--. 


C" -i 




CU K _ M _ 2 C 




■ C£7. 


n" . 



In (58) the case N > C has been presented. When C=N, Vj and in B'\ do not get C as multiplier in (58) and 
C appears as a multiplier of U j in A" 2 and A'\ in (56) and (57). The multiplier of matrices Uj and Vj 
concerning the services increases by one in each row block from third row block as the row number increases by 
one, up to the row C+l and it remains C in row blocks after that as given above. 



l 2 + cu x _ 2 



CI', 



The basic generator (59) which is concerned with only the arrival and service is Q B = A" 0 + A'\ + A" 2 . This 
is also block circulant. Using similar arguments given for Model (A) it can be seen that its probability vector is 

w'= (^<^'^< »^) where w is as seen in Model (A), where w= (^>i<Pi, §2<Pz, tyktVk*) an d me stability 

condition remains the same as in Model (A). Following the arguments given for Sub Cases (Al) and (A2) of 
Model (A), one can find the stationary probability vector for Sub Case (Bl) of Model (B) also in matrix 
geometric form. All performance measures in section 2.3 including the expectation of customers waiting for 
service and its variance for Sub Cases (Al) and (A2) of Model (A) are valid for Sub Case (Bl) of Model (B) 
with M is replaced by N. It can also be seen that when N = C the system admits Matrix Geometric solution as in 
Model (A). 

3.2.2Sub Case: (B2) C > N 

The infinitesimal generator Q B 3 2 °f me Sub Case (B2) of Model (B) has the same block partitioned structure 
given in (34) for Sub Case (A3) of Model (A) but the inner matrices are of different orders and elements. When 
C > N > M, the BMAP/M/C bulk queue admits a modified matrix geometric solution as follows. The chain X (t) 
describing this Sub Case (B2), can be defined as in the Sub Case (Bl). It has the infinitesimal generator Q B 3 2 of 
infinite order which can be presented in block partitioned form given below. When C > N, let C = m* N + n* 
where m* is positive integer and n* is nonnegative integer with 0 < n* < N-l. 
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0 


0 


0 


0 


0 
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0 
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In (60) the states of the matrices are listed lexicographically as 0, 1, 2, 3, .... n, .... Here the vector n is of type 
lxN^ifej andn = ((n, 0, 1, l),(n, 0, 1, 2)...(n, 0, 1, /q), (n, 0, 2, l),(n, 0, 2, 2)...(n, 0, 2, k 2 ),..., (n, 0, k*, 
l),(n, 0, k*, 2)...(n, 0, k* k kt ), (n, 1, 1, l),(n, 1, 1, 2)...(n, 1, 1, kj, (n, 1, 2, l),(n, 1, 2, 2)...(n, 1, 2, 
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fc 2 ),...,(n, 1, k* l),(n, 1, k*, 2)...(n, 1, k*, k kt ),..., (n, N-l, 1, l),(n, N-l, 1, 2)...(n, N-l, 1, /q), (n, N-l, 2, 
l),(n,N-l,2, 2)...(n,N-l,2,fc 2 ),...,(n,N-l,k* l),(n, N-l, k*, 2)...(n, N-l, k* k kt ) ) for n > 0. 
The matrices B"\, A"y for 1 < j < m* and A" l have negative diagonal elements, they are of order N £f=i fc; 
and their off diagonal elements are non- negative. The matrices A" 0 ,A" 2 j andA" 2 for 1 < j < m* have 
nonnegative elements and are of order N £?=i k t and the matrices A" 0 ,A" 1 andA" 2 are same as defined earlier 
for Sub Case (Bl) in equations (55), (56) and (57). Since C > N the number of servers in the system s equals the 
number of customers in the system L up to customer length becomes C= m* N + n*. Once number of customers 
L > C, the number of servers in the system remains C. When the number of customers becomes less than C, the 
number of servers again falls and equals the number of customers. Using the same matrices presented in model 
(A), for Q, Aj , Uj.Vj U, Q' and Q u in (6), (9), (10), (14) to (17) the partitioning matrices are defined below. 



matrix A" 2j 
-jNU N jNU N _, ' 
0 QN + 1)U N 
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jNU 2 
(jN + 1)U 3 



JNUi 
ON + V)U 2 



0 ••• 0N + N-2)U N QW + iV-2)f/ w _ 1 

0 ••■ 0 (JN + N - 1)U N . 

The matrix A 2mt is as follows given in (62) when C = m*N 
^(Nm*)U N (iVm ••• . ••• (Nm*)U 2 

0 (Nm*+1)U N ••• . ••• (JVm*+l)t/ 3 (/ 



A" 2 , 
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n* where 0 
' +1)U 2 



CU N 
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-- fi' - jU for 0 < j < C and Q 1C =Q as in Sub Cases (Al) and(A2). Then B"\, is defined as follows. 
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The matrix A" 11 
isC. 



n (65) when C = m*N + n* and 0 < n* < N-l. From row block n*+l, the multiplier of Uj 
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(65) 
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The basic generator for this model is also same as (59) which is concerned with only the arrival and service. 
Q B = A" 0 + A'\ + A" 2 . This is also block circulant. Using similar arguments given for Model (A) it can be 
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seen that its probability vector is w ' = (^<^<^< — ' where w = (cb^, 4 > 2 < P2: ^fc^fc*) ar, d the stability 
condition remains the same. Following the arguments given for Sub Case (A3) in section 2.2.2 of Model (A), 
one can find the stationary probability vector for Sub Case (B2) of Model (B) also in modified matrix geometric 
form. All the performance measures given in section 2.3 including the expectation of customers waiting for 
service and its variance for Sub Case (A3) are valid for Sub Case (B2) of Model (B) except M is replaced by N. 

IV. NUMERICAL ILLUSTRATION 

For the BMAP/M/C bulk models, the varying environment is considered to be governed by the Matrix Q 1 = 
| ^ J. Nine examples three for each are studied for the cases M = N =3; M = 3, N= 2 and M = 2, N = 3 
with the number of servers in each case as C = 2, 3 and 4. Matrix geometric results are seen for C = 2 and C = 3 
< M or N. Modified Matrix Geometric results are seen when C = 4 > M and N. 

The service time parameters of exponential distributions are respectively fixed in the two environments El and 
E2 as n x = 5 and \i 2 = .5 for single server respectively. 

For the case M=3, BMAP, the batch Markovian arrival process for El is given by Dq = ~ , = 



L. 32 



D\ -- 



F.12 .181 
L 08 .12V 



.121 



and BMAP for the environment E2 is given by Dq ■ 

«-E a. 



For the case M=2, Df for i=0 and i=l and D? for i = 0, 1, 2, 3 are as given above for the case M = 3 but it is 
assumed that Dl=[ 0 2 8 and D\ = [° °]. 

The bulk size service probabilities are given in table lfor the case when M = N = 3 for the two environment. 
For the case M =3, N=2 the probabilities of bulk service size 2 in El is fixed as .5 and of bulk size 3 in El is 
fixed as 0; and other probabilities are unchanged. 

Table 1: Service probabilities 

Environment 1 | P(size 1) | P(size 2) I P(size 3) I Environment 2 I P(size 1) I P(size 2) 



Service 



.5 



Service 



Thirty iterations are performed for all the models to iterate the rate Matrix R and the norms of convergence are 
recorded. Queue length probabilities and block size probabilities are calculated. Expected queue length and 
Standard deviation are presented. They show significant variations when M, N and C are changed. The 
probabilities of queue lengths and block sizes are presented in figures 1 and 2 for all the nine examples. 
Table2: Results Obtained For Six Matrix Geometric Models with Servers C=2, 3 and Three Modified 
Matrix Geometric Models with Servers C=4. 
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Figure 1: Probabilities of Queue lengths 



M=3=N M=3 / N=2M=2 / N=3 M=3=N M=3 / N=2M=2 / N=3 M=3=N M=3 ; N=2M=2 ; N=3 n(n>3)e 



Figure 2: Probabilities of Block Sizes 

V. CONCLUSION 

Two BMAP/M/C bulk service queues and their sub cases with randomly varying environments have been 
studied. The environment changes the batch Markovian arrival processes, the service rates, and the probabilities 
of bulk services. Matrix geometric ( modified matrix geometric) results have been obtained by suitably 
partitioning the infinitesimal generator by grouping of customers, environments, BMAP and PH phases together 
respectively when the number of servers is not greater than ( greater than) the maximum of the maximum arrival 
and maximum service sizes. The basic system generators of the queues are block circulant matrices which are 
explicitly presenting the stability condition in standard form. Numerical results for various bulk queue models 
are presented and discussed. Effects of variation of rates on expected queue length and on probabilities of queue 
lengths are exhibited. The decrease in arrival rates (so also increase in service rates) makes the convergence of R 
matrix faster which can be seen in the decrease of norm values. Bulk BMAP/PH/C queue with randomly 
varying environments causing changes in sizes of the PH phases may produce further results if studied since 
BMAP/PH/C queue is a most general form almost equivalent to G/G/C queue. 
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